Effect of an electric field on a floating lipid bilayer: a neutron reflectivity study 
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We present here a neutron reflectivity study of the influence of an alternative electric field on a 
supported phospholipid double bilayer. We report for the first time a reproducible increase of the 
fiuctuation amplitude leading to the complete unbinding of the fioating bilayer. Results are in good 
agreement with a semi-quantitative interpretation in terms of negative electrostatic surface tension. 

PACS numbers: 87.16.Dg, 87.15. Va,61.12.Ha 



I. INTRODUCTION 

In an aqueous environment, phospholipids self-assemble into bilayers; these 2D systems have been recently exten- 
sively studied, as they display complex peculiar behaviours P, and constitute the major component of the cell 
membrane [s*] . The strong influence of electric fields on lipid bilayers is interesting both from a fundamental point of 
view and for practical aspects, and has consequently received an important research effort. Large electric fields are 
also used to introduce macromolecules in vesicles or cells by formation of long-lived pores in the bilayer [1, Q , and also 
by endocytocis [6]. Smaller alternative electric fields are commonly used to destabilise lipid bilayers in order to form 
large unilamellar vesicles 0, Q in the so called electroformation technique. Recently Burgess et al. [9] have studied 
the effect of a static electric field (surface charge) on the structure of supported bilayers. 

Although it is a widely used technique, the mechanism involved in vesicle electroformation is not well understood 
yet. In a recent theoretical paper Sens and Isambert have suggested that a "negative surface tension" due to 
the electric field induces an undulation instability of the bilayer, and that hydrodynamical effects select the fastest- 
destabilised mode. These predictions are difficult to test in classical electroformation experiments. As a matter of 
fact, the initial system is usually a rather disordered lamellar phase, obtained by hydration of a dry lipid film and, 
consequently, the result is very polydisperse in both size and degree of lamellar order. Recently, Constantin et al. 
p"ll] have reported the first study of membrane destabilisation by an electric field on a more controlled system: a 
solid-supported fully-hydrated oriented lamellar phase (from 10 to 3000 bilayers). They observed for the first time 
individual peeling of the bilayers, but were not able to confirm or refute the existence of electrodynamical instability. 
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FIG. 1: Experimental set-up: the two identical silicon substrates (10 mm-thick) are in thermal and electrical contact with 
thermalisation blocks. The PTFE spacer ensures a 5 mm separation between the two blocks. 
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Because they are quite free of defects, double-bilayer systems are promising candidates to reach better understanding 
of the mechanism of vesicle formation in a well-defined geometry. Recent work has shown that they enable local study 
of the destabilisation of a single membrane [l^ . In this paper we focus on the first steps of the bilayer destabilisation 
by an electric field as detected by neutron refiectivity. 

II. MATERIALS AND METHODS 
A. Sample preparation 

Supported single and double-bilayers of DPPC and DSPC (1, 2-dipalmitoyl-sn-glycero-3-phosphocholinc and 1, 
2-distearoyl-sn-glycero-3-phosphocholine, Avanti Polar Lipids) molecules were prepared using the Langmuir-Blodgett 
and Langmuir-Schaefer techniques, as described in detail previously p^ . 

Bilayers were deposited on silicon substrates (highly doped 2-10 n.cm~^, 5x5x1 cm^^), polished to produce a high 
quality crystal surface with r.m.s. roughness < 0.2 nm (SILTRONIX, Archamps, F). After careful cleaning, the 
silicon blocks were made highly hydrophilic by a UV/ozone treatment as in (T3| . Lipid depositions were performed 
at a constant surface pressure of 40 mN/m, the Langmuir trough (NIMA, Warwick, UK) being thermalised at 20oC. 
Transfer ratios for the first layers were always high and of the same order as for samples described in [l^. A very 
precise adjustment of sample horizontality during the Schaefer dip, crucial for a good quality sample, led to transfer 
ratios always bigger than 0.96. After the deposition of the fourth layer, samples were constantly kept into water. 
They were closed while submerged in the trough, using a 5 mm-thick-PTFE spacer and a similar silicon substrate as 
a lid, and then placed into a water- regulated temperature chamber (see figure HI . 

Annealing was generally performed by progressively heating the sample to the fluid phase and then cooling it back to 
the gel phase. Then after careful caracterisation of the structure in the gel and the fluid phases, the adopted protocol 
consisted in applying to the supported fluid bilayer(s) an alternative electric fleld of fixed amplitude, decreasing 
gradually the frequency (from 200 Hz to 5 Hz). For each value of the frequency reflectivity data were collected. All 
measurements were made at 62°C to ensure sample fluidity. 

B. Neutron reflectivity: basic principles and measurements 

Specular reflectivity, R(q) , defined as the ratio between the specularly refiected and incoming intensities of a neutron 
beam, is measured as a function of the wave vector transfer, q = AttsiuB/X perpendicular to the refiecting surface, 
where 9 is the angle and A the wavelength of the incoming beam. R(q) is related to the scattering length density 
across the interface, p{z), by the relation [T4|: 

where |/5(g)| is the one-dimensional Fourier Transform of p{z). 

Specular reflectivity allows the determination of the structure of matter perpendicular to a surface or an interface. 
Experiments are performed in reflection at grazing incidence. Samples must be planar, very flat and with roughness as 
small as possible. Data were collected at the High Flux Reactor (HFR) of the Institut-Laue Langevin (ILL, Grenoble, 
F) both on the small angle diffractometer D16 (preliminary measurements) and on the high flux reflectometer D17, 
designed to take advantage of both Time-Of-Flight (TOF) and monochromatic methods of measuring reflectivity (dual 
mode instrument) [T5lll6l|. In this study measurements were all taken with the instrument in the TOF mode by using 
a spread of wavelengths between 2 and 20 A at two incoming angles (typically 0.7° and 4°) and with a resolution, 
defined by two choppers, of At/t = 1% and 5%, respectively, where At is the neutron pulse duration. The beam 
was defined in the horizontal direction by a set of two slits, one just before the sample and one before a vertically 
focusing guide. Neutron reflectivity from the same sample was measured at different temperatures monitored with a 
thermocouple (equilibration time 25 minutes, stability < 0.1°C, absolute precision < 0.3°C), in the water-regulated 
sample chamber already described in The useful q-range, before hitting the sample background, spanned from 
0.007 to 0.24 and was measured in about two hours, also after equilibration. The equilibration time varied 
depending on the step of temperature increase or decrease. For a step of 1°C it was of about 15 minutes. 
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FIG. 2: Reflectivity curves: closed circles (•) = double DSPC bilayer in CaCl2 before applying electric field and best flt with a 
double bilayer model (full line); open circles (o) = after applying a (10 V, 10 Hz) field and best fit with a single bilayer model 
(dashed line). In insert, refiectivity curves for the bare substrate before and after the experiment (and best fit with a silicon 
oxyde layer model). 

C. Analysis of reflectivity data 

Specular reflectivity is a technique able to resolve the stacking of several slabs each one with a certain thickness 
and a certain scattering length density. Data analysis requires a good knowledge of the sample and here it was done 
via model fitting. Born and Wolf gave a general solution, the so called optical matrix method to calculate the 
reflectivity from any number of parallel, homogeneous layers, which is particularly useful since any layered structure 
can be approximately described by dividing it into an arbitrary number of layers parallel to the interface, each having 
a uniform scattering length density. A reflectivity curve essentially reflects the sample density perpendicular to the 
substrate surface, or rather the square modulus of its Fourier transform. Since the phase is lost, data need to be 
fitted (e.g. with a slab model) to extract the density profile. In previous studies multiple contrast neutron 
measurements have determined within A precision the profile of adsorbed and floating bilayers, in the gel and the 
fluid phases. Each bilayer is resolved into outer-inner-outer slabs: outer slabs — heads (phosphocholine and glycerol 
groups), inner slab = tails (hydrocarbon chains). A water film of thickness D^, separates the adsorbed and floating 
bilayers. Taking into account the water layer between the adsorbed bilayer and the substrate and the silicon oxide 
layer, this leads to 9 slabs in total, each one having its own average thickness and scattering length density. Moreover, 
each interface between two slabs has a certain width, also called r.m.s. roughness. Since this enters into the density 
profile, it can be determined by fits of reflectivity curves, but it is not possible to determine whether it is a static 
intrinsic width or an average over the neutron beam's correlation length of temporal or spatial fluctuations. 



4 



8-1 0-' 




0.0 0.5 1.0 1.5 2 2.5 

q (nm-1) 

FIG. 3: Modification of tlie reflectivity curve witfi respect to the frequency of tlie electric field. A 5V alternative field is applied 
to a DSPC double bilayer in D2O; full circles (•) with f=100 Hz; stars (*) with f=50 Hz; open circles (o) with f=10 Hz. The 
lines are the best fits obtained with the values of the roughness given in the text. 

III. RESULTS 
A. Complete unbinding in a highly conductive solution 

The first result of this paper is the observation for the first time of the floating bilayer complete unbinding by 
application of an electric field. This effect was observed on two different samples under the same conditions: double 
bilayers of DSPC in a CaCl2 in D2O solution (0.35 mM.l~^, measured conductivity of 16.62 yitS.cm^^) at low frequency 
(10 Hz) and for a voltage amplitude higher than 5 V. Figure [2] shows reflectivity data for one of these double bilayer 
DSPC samples before and after applying the electric fleld. The curves correspond to the best fits obtained with a 
9-slab model and are in good agreement with previous multi-contrast neutron refiectivity experiments [T3| . 

After the unbinding we still observed the presence of the first bilayer, without any major change in structure, but 
slightly shifted away from the substrate (figure [2]) . After removing the bilayer by cleaning the substrate we could 
confirm that the silicon oxide layer was unmodified and therefore that there was no formation of porous silicon oxide 
(figure [5] in insert). The adsorbed bilayer might possibly be removed by application of a lower frequency or/and higher 
voltage field, but this could not be tested as it could have induced significant oxidation of the silicon. As a matter of 
fact, we had previously been able to determine that the formation of porous silicon oxide is very sudden and occurs 
when frequency becomes lower than 5 Hz in identical conditions (unpublished data). 

B. Progressive increase of fluctuations in a poorly conductive solution 

1. Double btlayer 

As this complete unbinding of the floating bilayer is very sudden, we tried to observe previous steps of the 
destabilisation in a low conductivity electrolyte such as D2O (xj^ — 0.34 /xS.cm"^). We performed experiments on 
a DSPC double bilayer in these conditions, and were able to observe a progressive modiflcation of the reflectivity 
curves (figure [3]). 
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• At high frequencies (f > 100 Hz), this modification appeared to be irreversible and involved mainly the first 
bilayer: the thin water layer (~0.3 nm, see table |l| initially separating the substrate and the first bilayer 
disappears (or becomes smaller than our experimental resolution), while the amount of solvent in the bilayer 
increases from 2% to 6%. For this first confined bilayer, assuming that the first water layer behaves as a 
dielectric, both effects could originate from the electrostatic pressure acting on the effective silicon oxyde-water 
layer-bilayer capacitance. Note that at the same time the second bilayer structure is unmodified. 

• At lower frequencies, modifications also involved the fioating bilayer roughness and were partially reversible. 
It is very difficult to determine which structural parameter (s) is/are responsible for the subtle modifications 
observed without additional information. Nevertheless, as previously mentioned no modification of the silicon 
oxide lay er was detected. The structure of the double bilayer can be well determined before applying the electric 
field [IjI, and the global structure of the bilayers (i.e. the slabs number and nature) seems to be unaffected by 
the electric field. It thus seemed physically reasonable to fit the refiectivity curves under electric field by varying 
mainly bilayers roughness and water layers thickness; average percentages of solvent within the bilayers were 
also allowed to vary, accounting for the induced formation of defects or pores in the membranes. 
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FIG. 4: Top bilayer roughness vs electric field frequency f at constant field amplitude lpp/2 = 5V: (•) DSPC double-bilayer in 
D2O (in insert, (■) same curve for a DSPC bilayer); (o) reversibility test after each frequency (no electric field) for the double 
bilayer. 

Within this framework, the dominant effect we observed was a progressive increase of the fioating (top) bilayer 
roughness, going from 0.6±1 to 1.3 ±1 nm. Again it was not possible to go to lower frequencies or higher voltages to 
avoid substrate deterioration, and in the frequency-range studied no complete unbinding occured in D2O. 

This roughness increase seemed to be partially reversible, as shown on figure ID Reversibility was systematically 
checked in the low conductivity solution (D2O); for beam time limitation reasons we were not able to check it in 
the CaCl2 conductive solution. A small increase of the thickness of water between the bilayers Dw,2, in the limit of 
the resolution of the experiment, was also observed. Finally, a progressive increase of the percentage of solvent in 
the first, substrate-bound bilayer was also observed, while the percentage of solvent in the fioating bilayer remained 
constant and very low. Whether this corresponds to the formation of pores - made easier by the higher tension of 
the supported membrane, or of permanent holes due to electrochemical phenomena close to the substrate is unsure. 
However the latter seems more likely as this increase of the solvent content was apparently not reversible. 
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The irreversible effect on the first bilayer at frequencies > 100 Hz, as weh as the reversible effect at low frequencies, 
well described by an increase of the top bilayer roughness, were observed on different DSPC samples. Nevertheless, 
as we were not able to check carefully the stability of the oxide layer in these cases we do not report corresponding 
data in the present paper. 



2. Single bilayer 

In order to check the effect of the electric field on the first, substrate-bound bilayer, we performed experiments 
on a single supported DSPC bilayer. Results are shown in figure H] (insert). They showed a progressive increase 
of the bilayer roughness, along with an increase of membrane-substrate distance. However no modification of 
the bilayer structure was seen, and the percentage of solvent remained very low (^--^2%), in agreement with a 
good-quality sample. We can thus think that, in the case of the double-bilayer, the solvent increase in the 
first bilayer is correlated to the fact that it is confined close to the substrate by the second, fluctuating bilayer. 
This seems consistent with the fact that when the second bilayer unbinds, the first one relaxes away from the substrate. 



IV. DISCUSSION 



The results reported above clearly show that an electric field can induce a complete destabilisation of the top bilayer 
of a double-bilayer at low frequency (~ 5-10 Hz). This very abrupt destabilisation follows a reversible increase of the 
bilayer roughness. It is tempting to describe this roughness in terms of thermal fluctuations, and wc intend to do so 
in this section. 



A. Bilayer fluctuations 



We consider the small thermal fluctuations of a single bilayer close to a substrate. Using common notations (see 
[ist for example), we note r = {x,y) the in plane coordinates and z the coordinate normal to the substrate and we 
call u{r) — u{x, y) the membrane position. Classical Helfrich's Hamiltonian can then be written at first order 
expansion in membrane's curvature as : 



n[u{r)] = / 



s 



U {u (r)) + (Am {r) f + 7 (Vu (r)) ' 



(2) 



where k is the bending modulus of the bilayer, 7 its surface tension and U the membrane-substrate interaction 
potential. Using a quadratic approximation for U {u{r)) it is possible to apply the equipartition theorem for small 
fluctuations: 

^'"^''^^^" [/" + 7"/+Kg^ 

Note that this expression is only valid when the fluctuation amplitude is smaller than the mean distance to the 
substrate, which seems to be the case in our experiments. For large fluctuations one would have to take into account 
the reduction of entropy caused by the substrate, leading to the so-called Helfrich repulsion [3. [Tsj. 

In their recent paper. Sens and Isambert ^1^ describe the case of a free membrane far from a substrate [U" = 0) 
and take into account the effect of hydrodynamics on the bilayer destabilisation. In particular, they showed that 
energy dissipation by solvent flow (viscosity 77) and by monolayer-monolayer friction (friction coefficient bfr) lead to 
the selection of a fastest-destabilised mode of wavelength q*: 

This length scale is in good agreement with the caracteristic size of vesicles prepared by electroformation. 

In our experimental study, the two bilayers are not free but interact with the substrate and the other bilayer. 
Bilayer equilibrium position without electric field is thus give n by a competition between hydratation repulsive forces 
and van der Waals attractive ones, as described in ref.f^. Tisj. 
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Using Eq. [3l it is possible to calculate the rms amplitude of fluctuations a^: 



^2 _ 



log — (5) 



27r\/A qi 

with A = 72 - AU"k, ql = 2^ and ql = 2±^. 

At the linear order, the effect of an electric field on the membrane can be expressed as a negative surface tension 
—Tel (with Tel > 0). Following the notations of ref. Tei can be written as: 

Tel = em£^m<5 (6) 

where em is the membrane dielectric constant, Em is the local electric field and 6 the bilayer thickness. This negative 
surface tension induces enhanced thermal fluctuations that can be described by replacing 7 by 7 = 7 — Tei in equation 
[3l Figure [5] shows the variation of ct as a function of the global surface tension 7 — Tef, here we have taken k ~ 40fcsT 
and U" ~ 10^^ J.m~^, which are classical values for a DSPC bilayer in the fluid phase [20j . 




Y-Tei (N.m-1) 

FIG. 5: Membrane rms fluctuation amplitude a as a function of global surface tension 7 — Tei, deduced from equation [5l (o): 
values probed during the experiments. 



Membrane fluctuations undergo a sharp divergence for a finite negative value of Tc = —\/U"k which is the signature 
of bilayer destabilisation. This unbinding occurs in a very narrow surface tension range, which could explain the abrupt 
destabilisation observed in our experiments. 
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B. Discussion in terms of electrostatic surface tension 



Let us now try to evaluate the electric surface tension equivalent to the field we have applied. It is possible to 
derive values of T^i usin g eq uation [51 if we know a, 7, k and U" . We measured a experimentally. Previous X-ray 
off-specular experiments [20| have shown on very similar samples (a fluid DSPC bilayer deposited on an hybrid OTS- 
DSPC bilayer) that k ~ 40 ksT and 7 ~ 0.5 mN.m"^. Using these experimental results we can fix U" to obtain 
the correct experimental value of a measured without any electric field. Following this method we obtain U" ~ 10^'^ 
J.m"''. This value is in good agreement with theoretical estimations using a van der Waals potential 0. It is then 
possible by using equation [5] to derive numerical values of Tei for the different fields applied, which are reported on 
figure [3 

Following ref. [13] , the electrostatic surface tension Tei (Eq. [6|) can be related to the electric field E applied to the 
whole cell by: 

r^^^e^^^^^EH (7) 

where Xm^ a-nd respectively the membrane and solvent conductivity. With ~ 2eo, — 0.34 x 10^"* S.m^^ 

for D2O (measured experimental value) and Tei ^ 10^^ N.m^^ it leads to Xm^ ~ 3 x 10"^*^ S.m~^ (equivalent to a 
membrane resistance R,„ ^ 0.2 Mfi.cm^), in good agreement with the value obtained by Purrucker et al. [2l|. In the 
high conductivity solution (CaCl2, xj^ = 16 x 10~* S.m^^), even if we found no value of the membrane conductivity 
reported in literature, it seems clear that we are well above the unbinding surface tension 7^. 

Field frequency seems to be an important parameter in membrane destabilisation. We were not able to detect any 
modification (reversible or irreversible) of the bilayer above 100 Hz whereas the effect becomes important below 10 
Hz (Fig. [Hand [7]). This is in good agreement with previous observations by Constantin et al. [ll|. It is also precisely 
the frequency range used in electroformation techniques. In order to understand more quantitatively this effect, let 
us try to estimate the local electric field acting on the bilayer. We can do so by using a simple equivalent circuit to 
represent the experimental set-up (figure [6]). 
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FIG. 6: A simple equivalent circuit for the system. 

The bilayer impedance Z,n is modelled by a resistance R,„ and a capacitance Cm in parallel (2]| . The substrate 
is described by a resistance Rsi (^ 1^^) and a capacitance attributed to the silicon oxide Csi02 (— 30 for a 1 
nm-thick oxide) in series. The electrolyte solution has a resistance Rs (— 92 kil for D2O). It is then easy to derive 
the negative electrostatic surface tension acting on the bilayer: 

r =ie (^VeI 

{r J S (l^ U^TmTif + LO^ [Tm + (1 -f n/3))^ 

where n is the number of bilayer(s), R the total resistance R — Rs + Rsi, P — and we have introduced the two 
characteristic times r„i = RmCm and — RCsi02- 

Using Rs = 92 kil and Csi02 — 30 /xF as fixed values, the best fit to the data gives R^ = 48 ± 1 kfi (~ 0.78 ± 0.2 
Mf^.cm^) and C™ = 0.1 ± 0.01 /iF (~ 0.01 ± 0.005 /iF.cm-^). The corresponding result is plotted figure[71 The value 
of Rm is in very good agreement with typical values of bilayer membranes resistance, around 0.9-1 Mfi.cm^ for a high 
coverage of the substrate (see for example Ref. [2l|, H^). However the value of Cm is surprisingly lower than values 
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Frequency f(Hz) 

FIG. 7: Tel as a function of frequency: (•) values derived from experimental data; dashed line = best fit using equation |8l full 
line — equation |8] with Km = 60 kfl and Cm = 2/iF. We have represented by a dashed line the frequency limit (f ~ 5 Hz) were 
silicon oxidation is observed. 

usually found for supported bilayers (~ 0.5 /iF.cm^^ j2ll. [2^). We have no explanation for this disagreement. Many 
effects could account for it, including the fact that the first membrane might be damaged by the field. Maybe our 
model is then too simple to give a good description of the system capacitive behaviour. We have drawn on figure [7| 
the curve obtained with a more consistent value of the membrane capacitance. In any case, we can notice that the 
frequency behaviour we observe is similar to previous observations by Constantin and coworkers 

V. CONCLUSIONS AND PERSPECTIVES 

In this work we have reported the effect of an alternative electric field on lipid double bilayers adsorbed on doped 
silicon substrates, as studied by neutron reflectivity. We clearly show that the electric field can lead to a complete 
unbinding of the floating bilayer in a high conductivity solution at low frequency. In a low conductivity solution, 
small reversible modifications of the reflectivity proflles are observed and described as changes of the rms roughness 
of the floating bilayer. These data could be interpreted with a fluctuation spectrum including a negative electrostatic 
surface tension term as suggested by ref. [10]. Electric field frequency seems to be an important parameter as no 
effect is observed above 100 Hz. This result is also consistent with previous observations [IJJ. We have tried to 
understand this effect with a simple electrokinctic model of the cell. 
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To the author's knowledge, this is the first time that the effect of an alternative electric field on the fluctuations 
of a membrane close to a wall is studied. We have clearly shown that it is possible to completly destabilise a single 
bilayer by using an electric field in a well defined geometry. This work should be useful in the future for a better 
understanding of phenomena such as vesicle electroformation. 

However, so far our results provide no information on the lateral length scale {q* in ref. [l3|) of the instability, 
which is one of the main predictions of the theoretical model [13] . We only measure the integrated root-mean-square 
amplitude a of fluctuations. In a near futur, we plan to perform off-specular reflectivity experiments, as described in 
ref. [20| . They could provide important information concerning the lateral length scale(s) of the changes involved, 
and will enable direct measurement of the fluctuation spectrum and the surface tension. 

Acknowledgments: The authors wish to thank Michel Bonnaud for technical assistance and the ILL for providing 
neutron beam-time. 



TABLE I: Structural parameter of a DSPC double bilayer in D2O before any electric field is applied: (sld) — scattering leng 
densities in 10~^ thicknesses and roughnesses in A. All values are obtained from the best fit (figure [2]). 





Silicon oxide 


Solvent 


Head Chains Head 


Solvent 


Head Chains Head 


thickness 

sld 
% solvent 
roughness 


11 ± 1 
3.4 ± 0.1 

1 ± 1 
3 ± 0.5 


3 ± 0.5 

100 
2 ± 1 


6 ± 1 36 ± 1 6 ± 1 
2.05 ± 0.1 -0.6 ± 0.05 2.05 ± 0.1 

19 ± 3 2 ± 1 19 ± 3 
3 ± 1 4 ± 1 3 ± 1 


19 ± 0.5 

100 

5 ± 1 


6 ± 1 36 ± 1 6 ± 1 
2.05 ± 0.1 -0.6 ± 0.05 2.05 ± 0.1 
19 ± 3 2 ± 1 19 ± 3 
3 ± 1 6 ± 1 3 ± 1 
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